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O 1 1 INTRODUCTION 

Compton scattering of cosmic microwave background pho- 
C^j tons off electrons in ionised intra-cluster gas leads to inten- 
IL* ' sity distortions along the line of sight (Sunyaev & Zel'dovich 
. \ 1972, 1980). The same scattering mechanism will also change 
■ the polarization state of the CMB if the radiation has a 
S_t ' quadrupolar anisotropy in the average rest frame of the 
electrons. Such a quadrupole can arise from several phys- 
ical effects, including: (i) transformation of a primordial 
quadrupole moment in the frame in which the CMB dipole 
vanishes (hereafter referred to as the CMB frame); (ii) aber- 
ration and Doppler effects in the transformation of the non- 
quadrupole moments of the primordial anisotropy due to 
electron bulk velocities; (iii) anisotropic scattering in the 
cluster; (iv) gravitational lensing from dynamic structures. 
The first three of these mech anisms were originally stud- 
ied by Sunyaev fe Zel'dov ich ( sunyaev & Zel'dovich 1980 



ABSTRACT 

We derive analytic expressions for the leading-order corrections to the polarization 
induced in the cosmic microwave background (CMB) due to scattering off hot electrons 
in galaxy clusters along the line of sight. For a thermal distribution of electrons with 
kinetic temperature k-&T c ~ lOkeV and bulk peculiar velocity V ~ 1000 kms -1 , the 
dominant corrections to the polarization induced by the primordial CMB quadrupole 
and the cluster peculiar velocity arise from electron thermal motion and are at the 
level of ~ 10 per cent in each case, near the peak of the polarization signal. When 
more sensitive measurements become feasible, these effects will be significant for the 
determination of transverse peculiar velocities, and the value of the CMB quadrupole 
at the cluster redshift, via the cluster polarization route. 
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general - polarization. 



Zel'dovich fc Sun y aev 198C ) , while the fourth was considered 



by Gibilisco Ql997p . 

Although the small levels of polarization expected from 
cluster effects fall below current experimental sensitivities, 
the prospect of detection with future CMB experiments, 
such as the Planck Surveyor satellite, has fuell ed a recent 
rise in interest in the topic. Audit & Simmons (1999) have 
(numerically) computed the frequency dependence of the 
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polarization arising from the transformation of the CMB 
monopol e for a cold cluster, and Itoh, Nozawa & Ko- 
hyama (1999) derived an analytic expression for the fre- 



quency dependence of this polarization mechanism, includ- 
ing the first two corrections arising from thermal effects. 
The claims from these two groups that the original results 
of Sunyaev & Zel'dovich w ere in error were met with a paper 
from Sazonov & Sunyaev ( 1999 ), who confirmed the earlier 
Sunyaev-Zel'dovich analysis, and derived analytic results for 
the frequency dependence of the polarization arising from 
the CMB monopole and primordial quadrup ole in the limit 
of a cold electron gas. Sazonov & Sunyaev (1999) also dis- 
cussed the polarization effect due to prior anisotropic scat- 
tering (which is second-order in the optical depth, r). 

In this paper we reconsider the question of polarization 
generation from clusters, in the limit where second-order op- 
tical depth and gravitational lensing effects can be ignored. 
We use a convenient, covariant representation of the linear 
polarization, which simplifies the geometric interpretation 
of our results, to derive simple analytic expressions for the 
frequency dependence of the polarization signal. Our expres- 
sions include leading-order thermal and kinematic correc- 
tions to the results derived by Sazonov & Sunyaev (1999), 



and confirm their results in the limit of low electron tempera- 
tures and bulk velocities. Our results for the thermal correc- 
tions to the polarization arising from the electron bulk veloc- 
ity differ significantly from those given by Itoh et al. (1999) 
We show that for typical cluster parameters, the leading- 
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order corrections to the polarization from the primordial 
quadrupole and the electron bulk velocity amount to reduc- 
tions in these effects of ~ f per cent at the peak of their 
signals. These corrections should thus be considered in any 
future analysis of forthcoming CMB polarization data in the 
directions of clusters, when such detections are permitted by 
improvements in receiver sensitivity. 

The paper is arranged as follows. In Section ^| we out- 
line our covariant procedure for calculating the scattering 
term, which appears in the the Boltzmann equation, for an 
unpolarized, but otherwise arbitrary, radiation field in the 
CMB frame. In Section we present our results for the 
polarization signal arising from the individual multipoles of 
the primordial CMB intensity, and discuss their relative im- 
portance for a typical cluster. We end in Section [| with our 
conclusions. We employ units with c = h = 1 throughout, 
and a H signature for the spacetime metric. 



2 POLARIZATION DUE TO ELECTRON 
SCATTERING OF AN ANISOTROPIC 
RADIATION FIELD 

We consider observations of a diffuse radiation field by ob- 
servers comoving with a timelike four- velocity field u a . We 
adopt such a Lagrangian (or I +3 covariant) approach since 
it provides arguably the most transparent description of the 
physics. For a recent review of the 1+3 covaria nt app roach 
in the context of relativistic cosmology see Ellis (1998). Rel- 
ative to «°, a photon with four- momentum p a is observed 



to have frequency v and propagation direction e a such that 
p a = v(u a + e a ). The photon propagation direction e a lies 
in the instantaneous rest space of the observer, so that 
u a e a — 0. We refer to geometric objects which are perpen- 
dicular to u a on every index, such cLS €■ ; ELS projected tensors. 
Since T homson scattering doe s not generate circular polar- 
ization ( Chandrasekhar 1950), we need only consider lin- 



early polarized radiation here. The linear polarization state 
of radiation with frequency v, which is propagating along 
e a , is most conveniently described by the symmetric linear 
polarization tensor V a b(v,e c ). The linear polarization is a 
trace- free, project ed ten sor which is also perpendicular to e a . 
Following Thorne (1980), we refer to projected tensors which 
are also orthogonal to e" as being transverse, and we denote 
the transverse, trace-free part of a rank-two tensor by the 
superscript tt, so that V a b = [V a b] TT ■ It is useful to intro- 
duce the (screen) projection tensor H a b = g a b — u a u b + e a e b , 
where g a b is the metric of spacetime. Hab is a transverse 
tensor which projects perpendicular to u a and e a . Making 
use of TLab, the properties of V a b{y, e c ) are summarised by 



Vab — V(ab), H a Vbc — Vab 



n ab Vab = o, 



(1) 



where round brackets denote symmetrisation on the en- 
closed indices. The linear polarization tensor is related to 
the Stokes parameters Q(v,e c ) and U(y,e c ) via 



Pn = 0/2, V12 = Vix = U/2, V22 = -0/2, 



(2) 



where the components of Vab are on a comoving orthonormal 
tetrad with the 3-direction aligned with e a . 

The generation of polarization through electron scat- 
tering is most simply described in the original rest-frame of 
the electron. Since we shall ultimately give our results in 



the frame in which the CMB dipole vanishes, we will need 
the transformation laws for the radiation variables under 
changes of the velocity field u a . If we consider an alternative 
four- velocity field u a , this can be written in terms of observ- 
able quantities in the u a frame as u a = 7(tt a + v a ). Here v a is 
a projected (with respect to u a ) four- vector which describes 
the relative velocity of an observer comoving with u a with 
respect to the u a frame. The quantity 7 is the Lorentz fac- 
tor associated with the relative motion, 7 = u a u a or equiva- 
lently 7~ 2 = 1 + v a v a . Similarly, we can express u a in terms 
of quantities in the u a frame, so that u a — j(u a + v a ), 
where v a is the relative velocity in the u a frame. Intro- 
ducing the projection tensor h a b = g a b — u a u b in the u a 
frame, and a similar projection in the u a frame, h a b, we 
have v a = —h b l v b /j. If we consider a given photon with 
four- momentum p a , the frequency v and propagation direc- 
tion e a observed in the u a frame are expressed in terms of 
the quantities in the u a frame via the Doppler and aberra- 
tion formulae: 

v = 71/(1 + e a v a ), (3) 
e a = [7(1 + e h v b )Y 1 (u a + e a ) - -y(u a + v a ). (4) 

For a given photon four-momentum p a , the screen projection 
tensor in the u a frame, TC a b, is related to TL a b by 

Hab = Hab - 2[i/(l + e c v c )]~ 1 p( a H h )dV d 

+ [i/(l + e c v c )]~ 2 p a pbHd 1 d 2 v dl v d2 , (5) 

while the linear polarization tensor in the u a frame, 
Vabi^,^), can be expressed in terms of the linear polar- 
ization in the u a frame as 



v~ z Vab{v, e) = v^H^H^Vd^ {v, e c ). 



(6) 



This transformation law ensures that the degree of polariza- 
tion V{v,e c ) = ^(2VabV ab )/I(v,e c ), where I(f,e c ) is the 
radiation intensity, is independent of frame for a given pho- 
ton momentum. 

The dynamics of the radiation is described by the colli- 
sional Boltzmann equation, which for the linear polarization 
takes the form 



C\v 3 Vab(v,e c )] = K a b(v,e c ), 



(7) 



where the Liouville operator acts along the photon path 
x a (X), p a {\) in phase space, with p a = da; a /dA, according 
to 



C[Vab{v,e c )} = H^HtWeVdM^^), 



(8) 



where V a is the spacetime covariant derivative, and 
K a b(v,e c ) describes the scattering. For CMB photons scat- 
tering off free electrons in hot clusters, the average photon 
energy in the rest frame of the scattering electron is small 
compared to the electron rest mass, so that electron re- 
coil can be neglected. Furthermore, the effects of induced 
scattering and Pauli blocking are negligible, so that the 
electron-radiation interaction can be approximated by clas- 
sical Thomson scattering. In this limit, the scattering term 
in the rest frame of a beam of electrons (which has four- 
velocity u a ) with number density n e , is given by (Challinor 
1999) 



b(i>,e c ) = n e ar(-V a b(i>, e c ) 
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(9) 



where or is the Thomson cross section. Here, I a b{v) and 
Sab(y) are projected (with respect to some general u a ) sym- 
metric, trace-free (PSTF) tensors which describe the inten- 
sity and electric pol arization quadr upoles respectively in 
a covariant manner flChallinor 1999| ). In equation (|), the 
tildes on I a b and £ a b denote that the intensity and polariza- 
tion quadrupoles are evaluated in the rest-frame of the elec- 
tron beam, and the transverse trace-free operations should 
be understood to refer to the u a frame also. For general £, 
we define PSTF intensity multipoles 



A, 



drie {Ae) I{v,e c ), 



(10) 



where A e = 4tt(-2)^!) 2 /(2^ + 1)! and dQ is the element of 
solid angle of e a , so that 



I(u, 



^2 l A, {v)e / 



(11) 



where Ae is shorthand for the index string a\ . . .ae and 
e A( = e ai ■ ■ ■ e ae . Angle brackets denote the PSTF part on 
the enclosed indices. The electric polarization multipoles are 
defined by 



2£(£- 



A« (£ + !)(£ + 2) 



dne {Af _ 2 V ai _ iae )(v,e c ), (12) 



for £ ^ 2. A complete multipolar description of the linear po- 
lariza tion requires the introduction of magnetic multipoles 
also ( Karnionkowski, Kosowsky fc Stebbins 1997 ; Seljak & 
Zaldari faga 1997| ; jChallinor 1999| ), but we shall not need 
their definition here. 

In this paper we shall only be concerned with the po- 
larization induced by the cluster in the limit of low optical 
depth. Furthermore, we shall ignore the effects of any pri- 
mordial polarization since these will be small compared to 
the effects of the intensity quadrupole. With these restric- 
tions, only the intensity quadrupole prior to the interaction 
with the cluster need be retained on the right-hand side 
of equation (|^). Transforming to the frame of an observer 
whose four- velocity is u a , the scattering term becomes 



K ab (u,e c ) = -^n c a T i> 2 Hi 1 H b l2 [I dl d 2 (v)] TT , 



(13) 



where the transverse trace-free operation refers to the u a 
frame, and we have used the transformation law for the scat- 
tering term: 



K a b(y,e) =Hi 1 Hb 2 K dl d 2 {v,e c ). 



(14) 



We shall compute the scattering term in the CMB frame, 
which has four- velocity it", so that the intensity dipole I a {v) 
vanishes. First we calculate the intensity quadrupole in the 
rest frame of a beam of electrons, in terms of the multipoles 
Ia c i v ) m the CMB frame. Making use of the invariance of 
I(v, e c )/v 3 , we have that 



IabW = T- 



15 



8tt ^ 





I 7 1 



[ 7 (1 + e<%)p+ 3 ) 



.(e c 



e )e( a e{,> j, 



(15) 



where v a is the relative velocity of u a with respect to u a , 
so that u a — 7(u" + v a ), and the integral is over solid an- 
gles of the projected directions e a in the u a frame. Next, 
we use this result for I a b(p) in equation ( |l3| ) to obtain the 
scattering term in the CMB frame due to a beam of elec- 
trons with number density dn e in its rest frame, and with 
velocity v a = -h%v b /-y relative to the CMB. Finally, wc ex- 
press dn e in terms of the number density of the beam in 
the CMB frame, using dn c = 7dn c , and integrate over the 
electron distribution. Assuming a non-degenerate, thermal 
distribution of electrons at temperature T e in the frame in 
which the electron bulk velocity vanishes, we have, in the 
CMB frame, 



dn c 



n c d 3 p c 



47rm3e c A' 2 (l/e e ) 



ex P [- 7 r(i + v a v a )/e e ], 



(16) 



where V is the projected, relative bulk velocity, F is the 
associated Lorentz factor, and d 3 p e is the element of 3- 
momentum. Here, C = feT c /m c is the dimensionless elec- 
tron temperature, m c is the electron mass, K?,(x) is a mod- 
ified Bessel function, and n c is the electron number density 
in the frame in which the bulk velocity vanishes. Note that 
the electron number density in the CMB frame is Fn e . 



3 POLARIZATION FROM THE MULTIPOLES 

In this section we compute the contribution to the scatter- 
ing term Kabiy, e c ) from each primordial multipole Ia c {v) 
in turn, expressing our results as expansions in the small 
parameters 6 C = fe B T e /(511 keV) and V = V~V a V a . 

In linear perturbation theory, where there is no non- 
linear scattering, there are no spectral distortions in the 
primordial radiation, so that 



I{y, e 



2v A 



pv/[k B T T (l+6 T )] 



(17) 



where T r is the all-sky average CMB temperature, and 
Sr{e c ) is the dimensionless temperature anisotropy. If we 
expand <?r(e c ) in PSTF multipoles TA e (Maartens, Ellis & 
Stoeger 1995] ), 



1=1 



(18) 



(recall there is no dipole in the CMB frame), then, to first- 
order in the primordial anisotropy, we have for £ > 1 



l(X)TA, 



where, i(x) = 



(19) 



with I(v) the intensity monopole (a Planck spectrum with 
temperature T r ), and x = f/(feT r ) the dimensionless fre- 
quency. In this manner, the frequency dependence of the in- 
tensity multipoles is contained in the scalar i(x) only, which 
proves to be very convenient for subsequent calculations. It 
should be noted that since 5r(e c ) has no monopole compo- 
nent by definition, no confusion should arise from our later 
use of the symbol r for the optical depth through the cluster. 
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3.1 The monopole contribution 

The contribution to the intensity quadrupole in the rest 
frame of the electron beam from the intensity monopole in 



the C1V D frame, 1(f), follows from equation 



15 



which evaluates to 

where /3 2 = —v a v a = -ti a !)„. Note that the contribution to 
Iab{v) from the CMB monopole is 0((3 2 ). To transform to 
the CMB frame, we require the result 



'^a 1 ^|; 2 [' ) ( C l« C 2>] TT = [7(1 +eV)] 2 [«{a«6)] 



(22) 



which follows from elementary manipulations. The trans- 
verse trace-free and PSTF operations on the left-hand side 
of equation (|2^) refer to the u a frame, while those on the 
right refer to the u a frame. It follows that the scattering 
term K a i,(v,e c ) due to a beam of electrons with relative ve- 
locity v a is proportional to lv( a v b < ) ] TT , so that the radiation 
is polarized along the normal to the plane containing v a and 
the line of sight e a . 

We use equations ( plj ) and (|2^) in the right-hand side 
of equation (|l^), and expand in powers of (3. The integral 
over the electron distribution is handled by an expansion 
in powers of V and Q a . The integrals are straightforward, 
but tedious, making them ideally suited to symbolic com- 
puter algebra packages (we use maple). The final result is 
an asymptotic expansion for the scattering term in the CMB 
frame due to the intensity monopole, of which the leading 
few terms are 



v K ab (is,e c ) = 



i MT [F (a F s) ]»^{^ 



Cv 



3F 



2(2F 2 +G 2 ) 



+ ~F(F* + 2G ) 



+ V> 



1 F +±(2F 2 +G 2 ] 



(23) 



where /i = — V a e a /V is the cosine of the angle between 
the bulk velocity of the cluster and the photon propagation 
direction. Here, F = a;coth(x/2) and G = xj sinh(x/2). The 
presence of the factor [V( a V{,)] TT ensures that the radiation 
is polarized along the normal to the plane containing the 
bulk velocity and the line of sight, as required by symmetry 
considerations. 

If we assume that the distance across the cluster is 
small compared to the characteristic length scale of the 
local spacetime curvature, then in the limit of low opti- 
cal depth the emergent polarization is simply V a b(v,e c ) = 
v 2 K ab (v, e c )t, where t is the light travel time across the clus- 
ter in the CMB frame. To obtain the degree of polarization 
V{v, e c ) from equation (p3[), one should multiply by t, take 
the ma gnitude (times V2), and divide by the emergent in- 



AI(i>,e c ) — n c OTtI(v) 



[@ e (F-A)+Vtj] 



Higher order terms i n the expansion of A7 



(24) 



have been 



1996 ; Challinor fc Lasenby 1999|; Etoh, Kohyama fc Nozawa 



1996 



I 



calcul ated recently ( Stebbins 1997 :_ Challinor & Lasenby 



1996 ; pSfozawa, Itoh & Kohyamal998; Sazonov & Sunyaev 



Identifying the effective optical depth through the 
cluster, r = n c aTt, in the limit of small C and V we recover 
the result for V (v, e c ) given as equation (22) in Sazonov & 
Sunyaev Q199SQ : 



• coth(a;/2)rV; 



^ eC ) = ^(e^T; 

where we have used 
2{V {a V b) ] TT [V^V b r T 



(25) 



(26) 

with V 2 = V (1 - m ) the square of the transverse velocity. 

The 0(6°, V°) and 0(&l, V°) terms in the curly braces 
on the right-hand side of equation (|3|) correct those given by 
Itoh et al. (1999). Itoh et al. also give an expression for the 
0(6 2 , V°) term; we include our result here for completeness: 

2 ri5 
A 2 

40 
7 



-F 



21(2F 2 + G 2 ) + ~^-F(F 2 + 2G 2 ) 
28 



F 4 + ^G 2 (11F 2 + 2G 2 ] 
F 4 + ic 2 (26F 2 + 17G 2 ) 



(27) 



In the Ray leigh- Jeans region, neglecting A/(z/, e c ) compared 
to I(y,e c ) w I(v), we find the degree of polarization 

1 



V(v,e c 



-rV 2 (l 

10 ' V 



69 c 



2v» + ±fe 2 e 



(28) 



which, as noted by Sazonov & S unyae v (199S), recovers the 
result of Sunyaev & Zel'dovich (1980) at lowest order. The 
0(Oc, V 1 ) term in the braces on the right of equation j23| ) 
is responsible for the weak dependence of the degree of po- 
larization on the absolute value of the bulk velocity, for 
a give n tra nsverse velocity Vt, as found by Audit & Sim- 
mons fll99S| ) through numerical calculations in the zero tem- 
perature limit. 

In Fig. [l] we plot the frequency dependence of the 
magnitude of the polarization, Hv)V(v,e c ), in units of 
2(kBT T ) 3 / (he) 2 , for the cluster parameters ksT c — lOkeV, 
V t = 1000 kms" 1 , t = 0.01. The direction of the bulk motion 
is chosen to maximise the 0(Q®, V 1 ) term in equation (p3|), 
so that fi = l/v3- For these parameters, the leading order 
correction to the polarization from the CMB monopole is 
thermal in origin, and tends to broaden the polarization sig- 
nal. In the Rayleigh- Jeans region the correction amounts to 
a reduction in the signal of ~ 10 per cent, which also fol- 
lows from equation (E8J) , while near the peak of the polarized 
signal (x = 4.85) the reduction is ~ 15 per cent. Near the 
peak, including the O(0 2 , V°) term in curly braces in equa- 
tion ( p3| ) gives a further 3 per cent correction upwards, while 
including the O(OcyV ) term (the leading order kinematic 
correction) gives a further correction upwards of only 1 per 
cent when the radial motion is towards the cluster. In the 



tensity I(v, e c ) + AI(v, e c ), where AI(u, e c ) is the change in 



intensity due to the scattering. Ignoring the small anisotropy 
in t he CMB frame, the leading order terms of AI(v, e c ) 
are ( |Sunyaev fc Zel'dovich 1972j ; Sunyaev fc Zel'dovich 198^ ) 



Wien region, where t he degree of polarization is large (Au 
Simmons 199S| ), thermal effects lead to a much larger 



dit 



fractional enhancement in the polarization. For example, in 
the Planck HFI channel at 545 GHz, the thermal enhance- 
ment of the polarization due to the monopole is nearly 30 
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Figure 1. Magnitude of the polarization [in units of 
2(A:B'jr) 3 /(fec) 2 ] due to the CMB monopole for a cluster with 
k B T c = lOkeV, V t = 1000 kms" 1 , r = 0.01, and fi = 1/V3- The 
dashed-dotted line is the contribution with only the O(@®,V ) 
term in curly braces in equation ( p^ ) included. The solid line 
also includes the O(Sl,V ) term, and the dashed line further 
includes the O(0 2 ,V°) term given as equation (|27j). The dotted 
line includes these first two thermal corrections, as well as the 
0(B°, V 1 ) term. 



per cent. However, for all of the Planck channels, the de- 
gree of polarization from the cluster parameters in Fig. [I] 
falls well below the predicted detector sensitivity. In the 
Planck HFI channels at 143, 217 and 545 GHz, the magni- 
tude of the polarization from the monopole as an equivalent 
fractional temperature change ST/T is ~ 0.014, 0.019, and 
0.070 uK/K respectively. These values should be compared 
to the predicted sensitivities of 3.7, 8.9 and 208 uK/K. Our 
conclusions concerning the feasibility of direct detection of 
cluster polarization with the Planck sat ellite are less opti- 
mistic than those in Audit & Simmons ( 1999 ). This is due 
to our use of rather more realistic cluster parameters, and 
the correct predicted equivalent thermodynamic tempera- 
ture sensitivities. 

The errors in our analytic expressions for the polariza- 
tion (and also the intensity change) decrease with decreasing 
cluster temperature and bulk velocity. For fixed parameters, 
the errors increase with increasing frequency. The latter be- 
haviour is due to the use of a Fokker-Planck expansion of the 
intensity with respect to the frequency change on scattering 
in the CMB frame, during the evaluation of I a b{v) [equa- 
tion (jl^)]. In practice this limits the range of frequencies for 
which the expansions are valid to x ^ 1//3 (recall that /? is 
the speed of the scattering electron relative to the CMB). 
For a typical cluster, the thermal electron motion is the lim- 
iting factor, giving a maximum dimensionless frequency x of 
the order of 10. Note that equation suggests a degree 
of polarization which behaves as x 2 for large x. However, 
this is not in conflict with the requirement that V(v,e c ) be 
bounded by unity given the limited range in x for which 
equation (B5|) is valid. 



3.2 The quadrupole contribution 

The contribution to the quadrupole in the rest frame of a 
beam of electrons from the primordial quadrupole in the 



Iab{v) = 



i[ 7 £(l + e%)] 
[ 7 (l + e d w d )]5 

~v ci )(e 



(29) 



Evaluating this integral is a little more involved than for the 
monopole calculation given in the previous subsection. An 
expansion in powers of /3 gives three distinct sets of terms. 
First, there are terms proportional to /i,^/i^r clC2 , where it 

will be recalled that h ab = g a b — u a Ub is the projection tensor 
into the instantaneous rest space of the u a frame. These 
terms transform to the CMB frame according to 

.TT - .TT 2 



[Tab 



1 



r cirr 

[V{a T b)c e I 
lV(aV b )] TT , (30) 



where a 



1 + e c v c . Second, there are terms proportional to 
c c , , which transform as 



'H C a 1 'Ht 2 lv{c 1 h'^ 1 2) v d2 T dld2 } 



\V(aT b ) c V 

•yet x ' 



+ 



(31) 



In equations © and @ the transverse trace-free and 
PSTF operations on the left-hand sides refer to the u a frame, 
while those on the right refer to the u a frame. Finally, there 
are terms proportional to Vi a v b \V cl v C2 t c1C2 , whose transfor- 
mation law follows from equation (|2^). It follows that, in 
general, the polarization due to the primordial quadrupole 
will be a superposition of six geometric objects: the five 
terms on the right-hand sides of equations ( |30[ ) and ([3l"|), 
with v a replaced by the cluster peculiar velocity V a , and 
a term V( a Vi,)V cl V C2 T clC2 . If the peculiar velocity of the 
cluster vanishes, the polarization tensor is proportional to 
[T a b] TT o nly. T he latter case was discussed by Sazonov & 
Sunyaev ( 1999] ) in the low temperature limit; it is straight- 
forward to show that the dependence of the degree of po- 
larization on the position of the cluster on the sky, for 
Vabiy-, e c ) oc [r a i,1 T T , is equivalent to that derived in Sazonov 



& Sunyaev (1999) by non-geometric means. 

The leading few terms for the scattering term in the 
CMB frame due to the primordial quadrupole are 

v 2 K ah {v, e c ) = ^-n c aTl{v)-^- 
1U e x — 



Y(M TT {l + e c [-6F 



+ i(2F 2 + G 2 )] + Vfi(-2 +F)X 



2[V( a T b)c e c ] 1 



(32) 



For completeness, we also give the 0(Q 2 ) term. This only 
has a [r a i,] TT geometric dependence, and so would appear in 
the curly braces in equation ([j2|): 

6 2 |-15F + ^(2F 2 + G 2 ) - yF(F 2 + 2G 2 ) 

}• (33) 



4 

+ 7 



F 4 + -G 2 (11F 2 + 2G 2 



by 



In the limit of small Q c and V, the polarization is given 



Vab(y,e) = jyTi(x)[T ab ] TT = — r[I ab (u) 



(34) 
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ues Ai = — 25 uK, A 2 = 23 uK and A3 = 2 uK used by 
Sazonov & Sunyaev (1999). The line of sight is chosen to 
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Figure 2. Magnitude of the dominant [T a (,] TT polarization [in 
units of 2(k^T I ) i / (he) 2 ] due to the CMB quadrupole for a clus- 
ter with k B T c = lOkeV, Vfi = 1000 krns" 1 and r = 0.01. We use 
the best-fitting COBE results for the quadrupole: Ai = — 25 |^K, 
A2 = 23 nK and A3 = 2 yiK. The dashed-dotted line is the con- 
tribution with only the O(0[?, V°) term in curly braces in equa- 
tion @ included. The solid line also includes the O(&l,V ) 
term, and the dashed line further includes the O(0 2 ,V°) term 
given as equation ^^j). The dotted line includes these first two 
thermal corrections, as well as the 0(0^, V 1 ) term. 



which also follows directly from equation Q13J) on setting 
v a = 0. This very simple result is equivalent to that derived 
by Sazonov & Sunyaev ( 199S| ), but the presentation here 



benefits from the fully geometric approach. We can write 
T a b m the form f r T ab = Y^i=i ^i n ° n b j wnere {V} are 
the (observed) eigenvalues of the temperature quadrupole, 
and {ni*'} are the corresponding unit projected eigenvec- 
tors. (For the metric signature employed here, T r T a h n^ = 
— AjTii 1 '.) For V a b{v,e c ) oc [r ab ] TT , it follows that the de- 
gree of polarization is maximal along the eigendirection with 
intermediate eigenvalue (say A3). Perpendicular to this di- 
rection, the polarization vanishes for directions such that 
(e a ria) 2 — (Ai — As)/(Ai — A2). For clusters at these posi- 
tions on the sky, the polarization signal due to the primordial 
quadrupole will vanish except for small corrections arising 
from the cluster peculiar velocity. The first such correction is 
given in equation (|32|). Including only this correction, the de- 
gree of polarization due to the primordial quadrupole along 
directions for which [r a b] TT vanishes evaluates to 

V = ^rK2[(A 3 - Ar)(A 2 - A 3 )] 1/2 , (35) 

which we have expressed as an equivalent temperature. In 
the Rayleigh- Jeans region, we find 



M TT [1-60. 



183 e2 
7 e 



2\V( a T b)c e c ] TT 



(36) 



Note that the 0(V) correction to the [r a f,] TT polarization 
term vanishes in this limit. 

In Fig. ^| we plot the magnitude of the dominant po- 
larization signal (that with [T a b] TT geometric dependence) 
due to the CMB quadrupole, for the cluster parameters 
k B T c = lOkeV, Vfj, = 1000 kms" 1 and r = 0.01. For 
the CMB quadrupole we take the COBE best-fitting val- 



lie along ni 3 ' to maximise the polarization signal. As with 
the polarization due to the CMB monopole, we see that 
the dominant correction to the polarization sourced by the 
primordial quadrupole is due to thermal effects, and these 
tend to broaden the signal. In the Rayleigh Jeans region the 
size of the fractional correction is the same as for the signal 
due to the monopole ~ 10 per cent [compare equations ( |28[ ) 
and (p6[)]. Near the signal peak (x — 3.83) the leading ther- 
mal correction leads to a reduction of ~ 14 per cent, but 
this is reduced to ~ 12 per cent when the next order thermal 
correction and the 0(V) kinematic term are included. The 
corrections are most significant in the Wien region; in the 
545 GHz Planck HFI channel the thermal terms increase the 
signal by ~ 60 percent. However, since the monopole contri- 
bution to the degree of polarization rises more rapidly with 
frequency (as x 2 ) than the quadrupole (as x) for large x, the 
quadrupole contribution is likely to be less significant than 
the monopole in the Wien regio n, for typical cluster parame- 
ters ( Sazonov fc Sunyaev 199S ). In the Planck HFI channels 
at 143, 217 and 545 GHz, the magnitude of the polarization 
from the quadrupole as an equivalent fractional temperature 
change 8T/T is ~ 0.016, 0.016, and 0.029 uK/K respectively, 
which, like the monopole contribution, is well below the pre- 
dicted sensitivities. 



3.3 The contribution of higher multipoles 

The £-th multipole of the CMB anisotropy gives rise to a 
quadrupole of 0(f3 ~ 2 ) in the rest frame of the scattering 
electron. It follows that the leading terms in the signal due 
to the CMB octupole (£ = 3) will be 0(V,O c ). Unlike the 
monopole, the octupole can give a thermal polarization sig- 
nal even if the cluster peculiar velocity vanishes (V — 0), 
due to the preferred directions in space that it defines. Since 
the octupole is of similar magnitude to the quadru pole [the 
COBE best-fitting valu es give r a t c |/|r a b ~ 2.3 (Btoeger, 
Araujo & Gcbbie 1997)], and given that the O(0 e ) term 



in the quadrupole contribution to the polarization can be 
a sizeable fraction of the total quadrupole signal, the oc- 
tupole contribution should potentially be included in future 
attempts to determine the primordial quadrupole at the lo- 
cation of distant (hot) clusters via the cluster polarization 
route ( Kamionkowski & Loeb 1997). 



The calculation of the polarization due to the CMB 
octupole follows that of the monopole and quadrupole. For 
the leading two terms we find 

v 2 K ab {v,e c ) = i neCTT / W _^l-^{[ Ta6o ^]TT (5 _ F) 
+ e c [r abc e c ] TT \&F - i(2F 2 + G 2 )] }.(37) 



For the hexadecapole (I = 4) the quadrupole in the rest 
frame of the electron is second-order in /3, so we might an- 
ticipate an O(0 C ) polarization signal. However, to this order 
the polarization from a beam of electrons with relative ve- 
locity v a is proportional to [T a bcdV c v d ] TT , which results in 
a polarization signal proportional to [r a b c dV c V d ] TT on inte- 
grating over the electron distribution. For reasonable values 
of the cluster peculiar velocity, we infer that the effects of 
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the CMB multipoles with I > 3 can be safely ignored, even 
for rather hot clusters. 



4 CONCLUSION 

We have presented a direct geometric method for determin- 
ing the polarization due to the scattering of anisotropic ra- 
diation off a thermal distribution of electrons, with non- 
vanishing bulk velocity, in the limit of low optical depth. 
Our results for the polarization due to the CMB monopole 
and quadru pole f ully agree with those obtained by Sazo nov 
& Sunyaev fll999| ) [see also Zel'dovich & Sunyaev ( |l98C| )] in 
the low temperature limit. Our results for the thermal cor- 
rectio ns to the monopole effect correct those given by Itoh 
et al. (1999). A careful inspection of their analysis suggests 



that the discrepancy arises from the omission of a rotation 
needed to relate the scattering basis in the rest frame of the 
initial electron and the CMB frame. We have shown that for 
typical cluster parameters the thermal corrections tend to 
dominate the kinematic corrections. Near the peak of their 
respective polarization signals, thermal effects tend to re- 
duce the magnitude of the polarization from the monopole 
and quadrupole by ~ 12 per cent. The corrections are much 
larger in the Wien part of the spectrum. We have also noted 
that the CMB octupole can lead to a polarization signal 
from hot clusters which is of similar magnitude to the ther- 
mal correction to the quadrupole signal. The effects reported 
here may be significant for future attempts to determine the 
CMB quadrupole at distant clusters and the tangential pe- 
culiar velocity, to an accuracy of a few per cent, via polar- 
ization observations. 
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